Measurement of the Mass Difference m(Df) — m(D + ) at CDF II 
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We present a measurement of the mass difference m(Dj~) — m(D + ), where both the Df and D + 
are reconstructed in the (pir + decay channel. This measurement uses 11.6 pb of data collected by 
CDF II using the new displaced-track trigger. The mass difference is found to be 

m(Df) - m(D + ) = 99.41 ± 0.38(stat.) ± 0.21(syst.) MeV/c 2 . 
PACS numbers: 13. 25. Ft, 14.40.Lb 



I. INTRODUCTION 



Meson masses are predicted by different models of quark interactions and the inter-quark potential. An- 
alytically, the spectrum of heavy-light mesons can be described in the QCD framework using the principles 
of Heavy Quark Symmetry and Heavy Quark Effective Theory P, • These theories state that in the limit 
of infinitely heavy quark mass, the properties of the meson are independent of the heavy quark flavor and 
that the heavy quark does not contribute to the orbital degrees of freedom. The theory predicts that up to 
corrections of order l/mh, c , m(B®) — m(B®) = m(Df) — m(D + )K$. Recently, lattice QCD calculations have 
also given their predictions for the meson mass spectrum 0, 0, @ ■ By measuring the masses of mesons pre- 
cisely, we narrow the range of parameters and approximations that theoretical models use to make predictions. 
For charm meson masses, a simultaneous fit 7] of all measurements including the mass difference between 
the and D + is used to compare experimental measurements with theoretical predictions. In this paper a 
measurement of the mass difference m(D+) — m(D + ) in the decay channels — > <f>ir + and D + — > </>7r + where 
4> — > K + K~ is presented 8]. The advantage of measuring the mass difference in a common final decay state 
is that many of the systematic uncertainties cancel. Gathering the large sample of charmed mesons used in 
this analysis is done using a novel displaced-track trigger, the Silicon Vertex Tracker (SVT) [jj, which enables 
recognition of the decay of long-lived particles early in the trigger system. 
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II. THE CDF II DETECTOR AND DATA SET 

The data used for this analysis were collected with the upgraded Collider Detector at Fermilab (CDF II) 
[Io| at the Tevatron pp collider. The integrated luminosity is 11.6 pb _1 at y/s = 1.96 TeV, taken during the 
period Feb - May 2002. These are the first physics-quality data from the Run II program. 

A. The CDF II Detector 

The CDF II detector is a major upgrade of the original CDF detector which last took data in 1996. The most 
important aspects of the upgraded detector for this analysis are the new tracking system and the displaced 
track trigger. CDF II, which is shown in Figure^ nas an integrated central tracking system immersed in a 1.4 
T solenoidal magnetic field for the measurement of charged-particle momenta. The innermost tracking device 
is a silicon strip vertex detector, which consists of three sub-detectors. A single-sided layer of silicon sensors, 
called Layer 00 (LOO) 0] is installed directly onto the beryllium vacuum beam pipe, at a radius of 1.7 cm. It is 
followed by five concentric layers of double-sided silicon sensors (SVXII) [13 located at radii between 2.5 and 
10.6 cm. The Intermediate Silicon Layers (ISL) |l3j | are the outermost silicon sub-detector systems, consisting 
of one double-sided layer at a radius of 22 cm in the central region, and two double-sided layers at radii 20 
and 28 cm in the forward regions. Surrounding the silicon detector is the Central Outer Tracker (COT) [l4j . 
a 3.1 m long cylindrical open-cell drift chamber covering radii from 40 to 137 cm. The COT is segmented 
into eight super-layers, each consisting of planes of 12 sense wires. The super-layers alternate between axial 
wires and wires with a ±2° stereo angle, providing three-dimensional tracking. This provides up to 96 position 
measurements on a track passing through all eight super-layers. A charged particle traversing the tracking 
volume deposits charge on nearby silicon micro-strips (clusters), and signals from the ionization trail in the 
COT are recorded by the sense wires (hits). Double-sided layers of silicon provide axial (r-ip) measurements 
of cluster positions on one side and z measurements via small-angle or 90-degree stereo information on the 
other. The LOO detector provides r-tp measurements only. COT information and SVXII r-ip information from 
the SVXII detector arc used in this analysis. 

B. Tracking Parameters 

CDF II uses a cylindrical coordinate system (r, ip, z) with the origin at the center of the detector and the 
z axis along the nominal direction of the proton beam. Tracks are fit to helical trajectories. The plane 
perpendicular to the beam is referred to as the "transverse plane" , and the transverse momentum of the track 
is referred to as px- In the transverse plane, the helix is parametrized with track curvature (C), impact 
parameter (do) and azimuthal angle p>Q. The projection of the track helix onto the transverse plane is a circle 
of radius r, and the absolute value of the track curvature is \C\ = l/(2r). The sign of the curvature matches 
the sign of the track charge. The do of a track is another signed variable; its absolute value corresponds to the 
distance of closest approach of the track to the beam line. The sign of do is taken to be that of p x d ■ 5, where 
p and d are unit vectors in the direction of the particle trajectory and the direction of the vector pointing from 
the primary interaction point to the point of closest approach to the beam, respectively. The angle ipo is the 
azimuthal angle of the particle trajectory at the point of closest approach to the beam. The two remaining 
parameters that uniquely define the helix in three dimensions are the cotangent of the angle 9 between the z 
axis and the momentum of the particle and zq, the position along the z axis at the point of closest approach 
to the beam. The two-dimensional decay length of a D meson L® y is defined as 

Y pD 

T D _ ±v • T m 
^ \P°\ (1) 

where is the transverse D momentum and X v is the vector pointing from the primary interaction vertex 
to the D meson decay vertex. We use the average beam position as an estimate of the primary interaction 
vertex. This is calculated for each data acquisition run. The transverse intensit y p rofile of the beam is roughly 
circular and can be approximated by a Gaussian distribution with a ~ 35 /jm 

mm 
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C. Trigger and Data Set 

CDF II has a three-level trigger system. The first two levels are implemented with custom electronics, while 
the third is a software trigger based on a version of the final reconstruction software optimized for speed. The 
most important feature of the trigger system for this analysis is its ability to recognize tracks and vertices 
displaced from the beam line. A brief description of this part of the trigger system follows. At Level 1 of the 
trigger, the COT provides information to the eXtremely Fast Tracker (XFT) that identifies tracks with 
Pt > 1-5 GeV/c. An event passes the Level 1 selection if the XFT finds a pair of tracks with opposite charge, 
such that each has px > 2.0 GeV/c, the scalar sum of transverse momenta pxi + PX2 > 5.5 GeV/c and angular 
difference A.tpe < 135°. The angle ip$ of a track is defined as the azimuthal angle of the track momentum as 
measured in super-layer 6 of the COT, which corresponds to a radius of 106 cm from the beam line. 

At Level 2, the SVT combines XFT track information with SVXII information. Tracks are refit using a 
linear algorithm, which provides improved (po and px measurements. The track impact parameter resolution 
is about 35 /im [15J, [l6( for tracks with px > 2 GeV/ c. An event passes Level 2 selection if there is a track 
pair reconstructed in the SVT such that each track has px > 2.0 GeV/c and 100 fim < \d \ < 1 mm. 

At Level 3, the full three-dimensional track fit using COT information is combined with SVT information. 
The Level 2 requirements are confirmed with the improved track measurements. The same tracks that passed 
the Level 1 selection have to pass the Level 2 and Level 3 requirements. In addition, it is required that 
the vertex of the two trigger tracks has L xy > 200 /im. The trigger requirements are optimized for selecting 
multi-body decays of long lived charm and bottom mesons. The optimization is done using an unbiased trigger 
sample to estimate the background rates and Monte Carlo simulated events to estimate the signal rates. 

Events gathered by the trigger system undergo final "offline" event reconstruction with the best available 
tracking algorithms. In the algorithm used for this measurement, the reconstruction begins with a COT 
measurement of the track helix. This version of the track is extrapolated into the silicon tracker, starting 
from the outermost layers and working inward. Based on the uncertainties of the track parameters, a road is 
formed around the extrapolated trajectory, and only silicon clusters found inside this road are added to the 
track. As clusters are added, the uncertainties on the track parameters are improved. For this analysis, only 
the r-ip information of the SVXII detector is used. 

III. MOMENTUM SCALE CALIBRATION 

The masses of the and D + mesons are measured from the momenta of their decay daughters, therefore 
it is crucial to calibrate the momentum measurements in the tracking volume. The main effects that are of 
concern in this analysis are a proper accounting of the energy loss in detector material and the calibration 
of the value of the magnetic field (B). Difficulties in accounting for energy loss in the tracking detectors 
come from an approximate model of the passive material. Uncertainties of the magnetic field are determined 
directly from the data. The momentum scale calibration for the tracking system is obtained by studying a 
sample of ~ 55,000 J/tj} — ► decays. An incorrect accounting for material in the detector description 

causes the reconstructed mass of the Jji\> meson to depend on the its px- Using an incorrect magnetic field 
value when converting track curvature into momentum causes the mass of the J /if) meson to be shifted. The 
calibration involves a two-step procedure. In the first step, the dependence of the J/ip mass on the transverse 
momentum is eliminated by adding material to the tracking volume description. After that, the magnetic field 
is calibrated by requiring that the reconstructed J/ij) — » mass be equal to the world average. 

A. Procedure 

The amount of passive material in the GEANT [l^ description of the CDF II silicon tracking volume 
is adjusted to eliminate the dependence of the invariant mass of the J/tp candidates on their transverse 
momentum, as demonstrated in Figure |21 The missing material is modeled with a layer of uniform thickness 
located just inside the inner shell of the COT; a layer of 0.56 ± 0.10 g/cm 2 eliminates the dependence of 
the J/tp — > mass on its px- This additional layer corresponds to roughly 20% of the total passive 

material in the silicon tracking system. Final state photon radiation causes a tail on the lower end of the J/i/j 
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mass distribution, which distorts (compared with a Gaussian distribution) the shape of the invariant mass 
distribution. The corresponding bias is calculated in bins of J/ip momenta and is taken into account when 
tuning the amount of passive material in the detector description. 

The magnetic field (B) is adjusted to bring the measured J /if) — ► fx" mass to the world average value of 
m(J/i/>) = 3096.87 MeV/c 2 0. The B field is calibrated to a value of 1.41348 ± 0.00027 T. The precision of 
the tuning procedure is limited by the number of J/tp decays available for the calibration. 

B. Tests and Cross-Checks 

Several tests and cross-checks are performed to verify the calibration. The J/ tjj invariant mass is checked for 
dependences on the z, ip and cot 6 coordinates of the decay in the detector. No significant residual dependence 
is found after the calibration is applied. The calibration method and parameters, the amount of missing 
passive material and the magnetic field value, are also cross-checked with other meson decays covering a range 
of invariant masses. As a check in the low momentum range, Kg — > n + ir~ decays are studied. The ir + 7r~ 
invariant mass distribution is presented in Figure |3| The Kg decays are also studied for dependencies on the 
radial position of the K s decay. No significant dependence is found for radii several centimeters inside the 
silicon detector. The mass of the K s is checked for run-to-run variations. No significant dependence on the 
run number is found. Cross-checks with high statistics, corresponding to several ten thousand signal events, 
are done with samples of D° — > K~ir + and D + — > K~ir + ir + decays presented in Figures 01 and 03 respectively. 
The D° decays are also checked for mass dependence on the pr of the D°. Since no particle identification is 
used, there is a reflection peak in the D° mass spectrum coming from the wrong assignment between kaon and 
pion hypotheses that can not be removed. The bias due to the reflection peak is estimated using a parametric 
simulation for every px bin separately and taken into account in Figure The ip(2S) — > J/i]jt{ + 'k~ decays 
are also reconstructed and the mass distribution is shown in Figure [7] Finally, a check in the region of higher 
momenta is done with T — > decays, presented in Figure |S] The reconstructed masses are compared to 

the world average values in Table HI We conclude that the calibration procedure described above accounts 
well for the energy loss in the silicon tracking volume, and applies to a range of reconstructed invariant 
masses. The calibration parameters quoted above are used when reconstructing the invariant mass of the 
and D + — ► <fm decays. 

One effect is found that is not completely corrected by the calibration. The distribution of the invariant 
mass of the J/ip as a function of the curvature difference between the two muons shows a slope, as seen in 
Figure [SJ This dependency indicates charge specific effects in the tracker, referred to as "false curvature". It 
also manifests itself in a difference in mass of the charge conjugates of the same meson. Misalignments in the 
COT, relative alignment of the COT to the silicon tracker, tilted wire planes and discrepancies between the 
COT axis and the magnetic field axis can cause such charge dependent false curvature effects. Parametrized 
corrections applied to track parameters improve the distribution shown in Figure |5J The charge asymmetry 
of the mass of charged mesons is not eliminated by these corrections. We do not correct for false curvature 
effects in the calibration procedure, but instead estimate the systematic uncertainty arising from the observed 
asymmetry. 

IV. D+ AND D+ SELECTION 

The Df and D + mesons are selected using offline reconstructed tracks through their decays to <t>iT + followed 
by the subsequent decay <fi — > A + A~ . To ensure good track quality, the tracks are required to have hits in > 20 
COT stereo layers, > 20 axial layers, > 3 silicon r-ip clusters and pr > 400 MeV/c. No particle identification 
is used in this analysis, and all mass assignments consistent with the assumed decay are attempted. 

The <f) candidates are selected by requiring two charged tracks, assumed to be kaons, which have opposite 
charge. The invariant mass of the track pair is required to be within 10 MeV/c 2 of the world average <f> mass. 
The detector resolution of the <fi mass is approximately 4 MeV/c 2 . A third track, assumed to be a pion, is 
added to the <fi candidate. To avoid using tracks from different interaction vertices, the separation along the 
beam line of all three tracks, the two kaon candidates and the pion candidate, is required to be < 4 cm. Any 
two of these three tracks satisfy trigger-like criteria using offline quantities: opposite charge, px > 2.0 GeV/c, 
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and 120 /im < \do\ < 1mm. The third track is only required to have |do| < 2 mm. No further requirements 
are placed on this track. 

All three tracks are constrained to a common vertex in 3 dimensions. To ensure quality of the vertices, 
the x 2 °f the vertex in the transverse plane satisfies x 2 ( r , f) < 7. The displaced track trigger preferentially 
accepts events with two-track vertices displaced from the primary interaction point by a few hundred microns. 
Adding a third track from the primary interaction pulls the three-track vertex toward the beamline, and the 
resulting L xy of the three track vertex is much smaller. To eliminate these background candidates, the L xy of 
the three track vertex is required to be larger than 500 /im. 

The helicity angle (Or) is defined as the angle between the 4> flight direction and the direction of the 
kaon momentum measured in the <f> rest frame. The <f> is polarized in this decay channel, so the helicity 
angle is expected to follow a cos 2 6h distribution for the signal, and a flat distribution for the background. 
Using sideband subtraction, we verify that the other selection requirements do not distort the shapes of these 
distributions, as demonstrated in Figure ITUl The helicity angle is required to satisfy |cos(##)| > 0.4. 

The requirements on the fit x 2 (r, tp), L xy and | cos(0h)\ have similar efficiencies. Individually, each require- 
ment is 90 — 95% efficient for the signal candidates, and rejects 40 — 50% of the background. It is unlikely 
to find two real Df/D + — > 07r + decays in the same event. If multiple candidates are found in an event, only 
the candidate with the highest | cos(9h)\ is considered. This procedure rejects another 9% of the underlying 
background. 

V. MASS FITTING AND SYSTEMATIC UNCERTAINTIES 

The invariant mass distribution of the K + K~ir + system is fit to two Gaussian distributions and a linear 
background. An unbinned maximum likelihood fit is used in which the widths of both Gaussian distributions, 
the mass of the D+ and the m(D+) — m(D + ) mass difference are allowed to float independently. Studies on 
both data and Monte Carlo simulation show that a linear dependence on mass is a good description of the 
background. Figure 1111 shows the likelihood fit superimposed onto the invariant mass spectrum. The \ 2 of 
the comparison of the likelihood fit to the measured mass spectrum is 127 for 118 degrees of freedom, and 
corresponds to a x 2 probability of 27%. The complete list of fit parameters can be found in Table ITT1 and the 
fit result yields: 

m(D+) - m(D+) = 99.41 ± 0.38 (stat) MeV/c 2 . (2) 

The two charmed mesons are produced either directly in the pp collision, or they are products of a B meson 
decay. The trigger preferentially selects mesons with large displacements of the decay vertex from the primary 
interaction point. Since the and D + mesons have different lifetimes, the fraction of directly produced 
J D + mesons to those coming from B meson decays is also different. Therefore, the momentum spectra 
of the two signals may differ, causing differences in the final state kinematics. This kinematic difference can 
produce a systematic shift in the measurement of the mass difference. Figure IT2"1 shows a comparison of the 
Pt distributions of the D+ (solid line) and the D + (dotted line). The spectra are very similar, and we expect 
small systematic uncertainties. 

A. Discussion of Systematic Uncertainties 

The systematic uncertainties are summarized in Table IIIII and will now be discussed in order of decreasing 
size. The largest single systematic uncertainty comes from fitting. To estimate the systematic uncertainties 
due to background modeling, the results of fits with different background models are compared. One model 
used in this comparison is a linear combination of orthogonal polynomials. Another model consists of two 
piecewise linear functions that meet at a point, which is varied between the D + and mass distributions. 
A systematic uncertainty of 0.08 MeV/c 2 on the mass difference is assigned based on the variation of the 
fit result when these different models are used. The systematic effect due to signal modeling is studied by 
excluding regions of the Df and D + signals from the fit. In this case, a fraction of the variation of the fit 
result is caused by changing the statistics of the sample used. This contribution is estimated by comparing 
statistical uncertainties of the fit result with regions excluded to that of the fit result with no modification. 
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After estimating the statistical contribution of the variation of the fit result, the systematic uncertainty due to 
signal modeling is estimated to be 0.12 MeV/c 2 . These two systematic uncertainties are added in quadrature 
and a systematic uncertainty of 0.14 MeV/c 2 due to fitting is obtained. This is the largest single systematic 
uncertainty. 

To estimate the systematic uncertainty introduced by sample selection requirements, the requirements on 
X 2 (t, y), L xy , cos,6h, and duplicate rejection are individually varied. Fit results were compared to estimate 
systematic effects for individual selection requirements. A fraction of the variation in the fit result is caused by 
statistical effects due to changing the sample composition when the selection requirements change. As before, 
the statistical contribution to the fit result variation is estimated from the change in the statistical uncertainty 
of the fit result. The only relevant selection requirement which exhibits a statistically significant effect is the 
cut on the x 2 (r, ip) variable. This variation of the mass difference is traced to an enhanced background around 
the D + mass for small values of the x 2 ( r ; f) variable. The effect is estimated to cause a systematic uncertainty 
of 0.11 MeV/c 2 . 

The systematic uncertainty due to the momentum scale determination is estimated by analyzing a kine- 
matically similar decay. A GEANT study is done to determine how the uncertainty on the mass difference 
measurement would scale with the absolute uncertainty on the D + — > Ktttt mass due to momentum scale 
variations, and shows that the uncertainty on the mass difference corresponds roughly to 11% of the absolute 
uncertainty on the D + mass. The world average mass of the D + meson m(D + ) — 1869.4 ± 0.5 MeV/c 2 is 
compared to our measurement of m(D + ) = 1868.65±0.07 MeV/c 2 obtained in a sample of D + — ► K~n + Tt + de- 
cays, using the same calibration procedure. To determine the absolute uncertainty of the momentum scale, the 
uncertainty of the world average (0.5 MeV/c 2 ), the statistical uncertainty of our measurement (0.07 MeV/c 2 ) 
and the difference between the two measurements (0.75 MeV/c 2 ) are added in quadrature. The sum in quadra- 
ture is then scaled by the factor obtained in the Monte Carlo study, and the systematic uncertainty of the 
momentum scale determination is estimated to be 0.10 MeV/c 2 . 

The mass difference is also sensitive to detector effects that are not corrected for by our calibration, namely 
false curvature effects. These effects are expected to cancel in the measurement of the mass difference. As 
explained in the calibration section, empirical corrections of the track curvature do not completely eliminate the 
asymmetry of charge conjugate states. By comparing fit results with and without these empirical corrections, 
the systematic effect of uncorrected tracking effects is estimated to be 0.06 MeV/c 2 . 

The accuracy of the momentum scale calibration is limited by the size of the J/ip sample. The systematic 
uncertainty on the mass difference from this limitation is estimated by individually varying the amount of 
material and the magnitude of the magnetic field by their statistical precisions. The two systematic effects 
are added in quadrature to obtain a systematic uncertainty of 0.03 MeV/c 2 due to the calibration procedure. 

Finally, an explicit check is done for a systematic uncertainty caused by the difference in px spectra of the 
Df and D + shown in Figure^l The events were re-weighted in the fit to make the spectra identical and the 
systematic effect on the mass difference is found to be negligible. 

The total systematic uncertainty of the measurement is estimated by combining the above systematic 
uncertainties in quadrature, and is found to be 0.21 MeV/c 2 . 

VI. SUMMARY 

The difference between the mass of the Df meson and the D + meson is measured using 11.6 pb _1 of data 
collected by CDF II and is found to be 

m(D+) - m{D+) = 99.41 ± 0.38(stat.) ± 0.21(syst.) MeV/c 2 . 

The result is in agreement with the current world average and the most recent Babar publication of 
(98.4 ± 0.1 ± 0.3) MeV/c 2 with a comparable uncertainty. 
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Decay Mass [MeV/c 2 ] PPG [MeV/c 2 ] 

K° s -^Tv + ir~ 497.36 ±0.04 497.672 ± 0.031 
T -^fJ, + n~ 9461 ±5 9460.30 ±0.26 

D° -^K~it + 1864.15 ±0.10 1864.5 ±0.5 
D + -^R--k + tv + 1868.65 ±0.07 1869.4 ±0.5 

ip{2S)-* J/tpTV+n- 3686.43 ±0.54 3685.96 ±0.09 



TABLE I: Table comparing measured masses of mesons reconstructed using the described calibration parameters and 
corresponding PDG averages. Uncertainties on reconstructed masses are statistical only. 



Parameter Value 

5m [MeV/c 2 ] 99.41 ± 0.38 

m{D s ) [MeV/c 2 ] 1968.4 ±0.3 

a(D s ) [MeV/c 2 ] 8.4 ±0.2 

a(D+) [MeV/c 2 ] 7.3 ±0.3 

f(D s ) 0.65 ±0.01 

/(£>++ D„) 0.37 ±0.01 

background slope [1/GeV/c 2 ] -7.3 ±0.7 

X 2 /NDF 126.7/118 (27.9 %) 



TABLE II: Table of likelihood fit parameter results corresponding to Figure [TT1 The x 2 , number of degrees of freedom 
(NDF) and corresponding probability are also listed. The parameters are the mass difference (5m), the mass of the 
Df meson, the mass resolutions (a(D a ), a(D + )), the fraction of signal events (f(D s ) , f(D + + D a )) and the slope of 
the background. 



Effect 


Syst.[MeV/c 2 ] 


fitting (signal ± background) 


0.14 


event selection 


0.11 


momentum scale 


0.10 


tracker effects 


0.06 


calibration procedure 


0.03 


Total 


0.21 



TABLE III: Table of systematic uncertainty estimates for the mass difference. The total uncertainty is the quadratic 
sum of the individual uncertainties. 
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END WALL 
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FIG. 1: Quadrant view of the CDF II integrated tracking system. The Central Outer Tracker (COT) and silicon 
subdetectors form an integrated tracking system. 
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FIG. 2: Dependence of the J/ip mass on the pr of the J/ip. The open squares show the mass dependence for tracks with 
no energy loss corrections. Open triangles show the result after applying the energy loss for the material accounted for 
in the GEANT description of the detector. Open circles account for the missing material modeled with the additional 
layer. Filled circles show the effect of the B field tuning in addition to accounting for all the missing material. 
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FIG. 3: Measured -k + tt mass distribution. A Gaussian distribution and a linear background are fitted to the mass 
spectrum. 
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FIG. 4: The K ir + mass distribution of the reconstructed D° candidates. A Gaussian distribution for the signal and 
a broad Gaussian distribution for the background are fitted to the mass spectrum. 
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FIG. 5: The K ti + tt + mass distribution of the reconstructed D + meson candidates. A Gaussian distribution for the 
signal and a broad Gaussian distribution for the background are fitted to the mass spectrum. 
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FIG. 6: The dependence of the D mass on its transverse momentum. The hollow points show mass values before 
any corrections are applied; the solid points show the dependence after the calibration (energy loss and B field). The 
systematic bias due to background modeling has been subtracted. 
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FIG. 7: Measured /i + fi -k + tt mass distribution for ip{2S) candidates reconstructed in the J/ip tv + -k decay. A 
Gaussian distribution and a linear background are fitted to the measured spectrum. 
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FIG. 8: Measured T — > fi + n mass distribution. Three Gaussian distributions and a linear background are fitted to 
the mass spectrum. 
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FIG. 9: Dependence of the J/ip mass on the difference of the absolute values of the curvature (C) of the positive and 
negative muon. This distribution shows a small charge dependent effect that are not corrected for in the calibration. 
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FIG. 10: Sideband subtracted and sideband distributions of the cosine of the helicity angle of the D + candidates (left) 
and Dt candidates (right). Candidates with \cos(6h)\ < 0.4 are rejected in the selection. 
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FIG. 11: Measured K K n mass distribution compared to the unbinned likelihood fit. 
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FIG. 12: Sideband subtracted distributions of the pr of the candidates (solid) and D + candidates (dots). Both 
distributions are normalized such that the sum of the bins add up to one. 



